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Abstract 

We report calculations which suggest that there is a physically-important 

four-atom exchange process in bcc He and thus an important four-spin term in 

the exchange Hamiltonian. A simple, mean-field analysis of this Hamiltonian 

3 

appears to account for a number of the perplexing properties of bcc He. An 
understanding of other properties may require treatment of the exact four-spin 
term. It is our hope to stimulate such effort by this Letter. 
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Wo report the results of calculations which suggest there is a physically 
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important four-atom exchange process in bcc solid Ho. The process gives rise 
to a four-spin term in the effective spin or exchange Hamiltonian with an 
exchange energy comparable to the nearest-neighbor two-spin term. A simple- 
minded mean-field treatment suggests this four-spin term could lead to a temp- 
erature-dependent exchange frequency which offers partial insight to the several 

3 1 

perplexing properties of bcc solid He. 

To facilitate discussion we define the exchange Hamiltonian including 

2 

pair, triple, and the important cyclic quadruple exchange: 

<i> - - 

H ex * " 2 CJ 1 - CJ 112 + 3J mi,23 3 X ± ' 
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Tlie first two-spin term involves nearest-neighbor spins (the "(l) 1 ’ over the sum), 
while the second involves next-nearest neighbor spins (the "(2)"). Finally the 
one four-spin term involves four atoms located at the corners of the rhombus, 
lying in the (110) plane, whose sides IJ, jk, k % H are first-neighbor dis- 
tances (the subscript "llll" in J 1]L11 23 > and whose diagonals ik, j 1 are second- 
and third-nearest neighbors (the 23 after the comma) . 

Later in this letter we will present a summary of the calculations for 
the various exchange constants. The relative large size of 23 appears to 

arise from the generally larger amplitude of atomic vibrations associated with 
transverse phonons. Moreover the ratio of 23^ J 1 is ^ urther enhanced by 

making the atomic wavefunction less localized than is indicated by most calculations. 
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of the ground state energy, which are relatively insensitive to the degree of 
localization, 

3 

First we want to briefly discuss the anomalies in He solid within the con- 
text of a mean-field approximation for (1) in which, for simplicity, we ignore 
the next-nearest-neighbor exchange terms except to indicate those situations 
where they may play a noticeable role. Then (1) reduces to 


( 1 ) _ 

H = - 2? (T) E I • I 
i<3 3 


where 


*!<*> = J, - 6J 112 + 3J 111123 (l + 4 <? • !>,). 


( 2 ) 


(3) 


In (3) the term (i • 1*)^ corresponds to the thermal expectation value of the 
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scalar product of nearest-neighbor spins. Since we expect He solid to become 
anti-ferromagnetic at low temperatures, then si • 1)^ will change from zero at 
high temperature to something like -(3/4) in the ordered phase. Taking a simple 
analytic form for this switching behavior, and using numbers we have calculated, 
we might approximate (3) above the transition. by 


^(T) w - 0.65 4 0.4 tanh (2/T) (4) 

where 3^ and T are measured in mK. (The factor 2 in the tanh (2/T) is an esti- 
mate of the transition temperature in the absence of J 23*^ 

3 

There are several puzzling anomalies in the He solid data which we discuss 
in terms of (£) and (4). We restrict ourselves to a single density, corresponding 
to a molar volume of 24.0 cc/mole. 

High temperature results (T > 20 mK) : In this regime, is a temperature- 

3 

independent constant. High temperature susceptibility measurements indicate 

it is negative (anti-ferromagnetic), with its best value of -0.65 mK coming 

4 

from the pressure measurements. When the pressure measurements are extended 

5 

to finite magnetic fields, it appears necessary to assume an additional 
ferromagnetic next-nearest neighbor interaction of strength 4, /J^ ~-0.2. 
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All of the high temperature c!nta can he reconciled within such n two-parnmetor 

^ 6 7 

model, in which 3^ is still approximately -O.GS mK, and is ferromagnetic. ' 

6 8 1 

Ratios in the vicinity of -0.2 have been theoretically calculated. * ’ In the 
present work we can only specify +0.03 > Jg/J^ > -0.2 because of 10 % or so 
uncertainties in the various exchange frequencies, and large cancellations in 

J 2 “ 4 J 112 + J llll 23 ^ h ° hieh tem P erature value 0:C J 2 (T) ) . 

Low temperature results (T<20 mK) : 

(i) Transition temperature . A constant first-neighbor exchange frequency of 
-0,65 mK corresponds (via T^=2.75|j|) to an anti-ferromagnetic transition 
temperature of 1.8 mK (a bit higher if 3g>0), whereas the first observed 

9 

transition in the solid is 1.1 mK. Clearly (4) would produce a lower transi- 
tion temperature. 

(ii) Specific heat . A high temperature expansion for the specific heat in 
powers of T _1 using an anti-ferromagnetic Heisenberg Hamiltonian results in a 

T term with a positive coefficient (even larger if 3" > 0) ; whereas the 
9 10 

experiments * indicate a negative one. A temperature-dependent exchange such 
as (4) can result in a negative coefficient, although our crude agrument appears 
to give too small a coefficient. 

Near the transition temperature of the solid the specific heat has a pro- 
9 10 

nounced bump ' near 2 mK prior to the entropy anomaly at 1.1 mK. Our mean 

field model cannot explain this structure. On the other hand, the original 

Hamiltonian (1) we propose may be able to explain it. Specifically we propose 

the phase transition theory community the following question: does H (1) 

cx 

with its four-spin term have the same behavior as a Heisenberg Hamiltonian in 

the vicinity of the (depressed) transition temperature? Or can it exhibit 

3 

structure consistent with that observed in bcc He? 

(iii) Spin Diffusion . As additional evidence that something extraordinary is 
occuring at about 2 mK in the solid, there is the fact that the spin diffusion 
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coefficient D drops sharply by over 30% for decreasing temperature about 2 niK,^ 
leading once again to the idea that the effective J («D) is strongly temperature 
dependent. 

(iv) phase boundary . In terms of a mean-field model, an external magnetic 
field would favor a positive (l • 7)^ in (3) and hence increase | | . Pre- 

sumably this effect would increase the magnetic field at which the spin-flop 

to paramagnetic-solid transition occurs over that predicted by a Heisenberg 

12 

Hamiltonian. Only preliminary experimental data exist for the H-T phase 
boundary. 

Finally we turn to the actual calculation of ^ 2 3 the ot ^ er exchange 

13 

frequencies. There is an extensive literature^ which underpins our calculation. 

The exchange frequencies can be written in terms of integrals involving the 

3 

Jastrow-Gaussian ansatz for the many-body wave function in solid He: 


♦ = TT e“ iA(r i " V TT f ( t )j 
1 i<j J 


(5) 


where r^ and specify the atom and lattice site positions, respectively. 

The scalar coefficient A is proportional to the (Einstein) frequency 

describing the harmonic motion of an atom about its lattice site. In our calcu- 

14 

lation we generalize A to a matrix in order to take into account a more 

15 

realistic description of the atomic motion based on a self-consistent calculation 

of the phonon frequencies. Accordingly the A . „ listed in Table I for the four 

eix 

largest exchange frequencies reflects phonon frequencies which predominate in 
that exchange process. 

Hie Jastrow function f(r .) is roughly a step function used to describe 

the short-range correlations, by preventing two atoms from approaching closer 

than their hard-core diameter. The energy of solid He can be written in terms 

16 

of integrals involving l|f. Various calculations which in effect minimize the 
energy with respect to functional form and parameters of (5), yield very similar 
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Jastrow functions. On the oti v hand the energy is not very sensitive to A. Never- 

g 

theless, it now appears on the basis of recent Monte Carlo calculations, that 
A may be 307o smaller than the range of values previously accepted. We exploit 
this possibility in our work. 

In terms of the wave function (5) the exchange frequency for a given ex- 
change p of atoms is 

2 

J -f) e -A d P /4 I (A,{#}). (6) 

pp p ’ 

where d^ is a distance in configuration space between the original and exchanged 

atoms. Values of d^ in units of the nearest neighbor distance H are given for 

various exchanges in Table I. The attempt frequency for the exchange process, 

within the Einstein approximation for the atomic motion, is given by 

0 = 0l/2M) d (A 3 /tt)“. 

P P 

Tiie exchange frequencies are dominated by the other two factors in (G) . 
Gaussian overlap . The exponential factor is a measure of the overlap of the 
wavefunctions of the original and exchanged atoms. Note that it is a very 
sensitive function of the frequencies of the exchanging atoms. In Table I A^^ 
for next-nearest-neighbor exchange is smaller than for nearest-neighbor which 
is a measure of the relative ease of distortions for exchange in (100) as 
opposed to (111) directions. 


Hard-core effects. The final factor I corrects for the fact that the Gaussian 

p 

overlap overestimates the probability for atoms to exchange by not excluding 
routes of exchange forbidden by the hard-core repulsion of the atoms. This 
effect is most important for pair exchange where straight-line exchange would 
have the atoms sitting on top of each other. In this case 1^ includes the 
dominant factor exp ^ C 2 ) jr where a is the hard-core radius. This strong 
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dependence on A does not occur in I for throe- and four-atom exchange, whore 

P 

tho short-range correlation effects are largely geometric in nature. 1 The 
values for in Table I are extrapolated from the results of Ref. 1 based on 
Monte Carlo integrations (for the Einstein case) , with the exception of the 
Inst (and crucial) entry which has been estimated on the basis of the experience 
gained from three- atom exchange. 

In the next to the last column of Table I are the results for J . It is 

P 

clear that J is too small but so is J, by at least a factor of 20, if it 

is to agree with the high temperature measurements. At this stage we take 

8 

advantage of the fact that the more accurate Monte Carlo work suggests that 

the true A values tr.ay be 30% or so smaller than those consistent with Ref. 15 

on which our work is based. Accordingly we decrease A, or more accurately, 

rescale downward che phonon frequencies until the value of + 3J 1111 23^ 

[see (1)3 agrees with the exchange frequency deduced from high temperature 
17 

measurements. The resulting numbers are shown in the last column. We note 
that 9 1 J mi 23 I is ™ K > the factor used in (4). As a result of these cal- 

culations we argue that this four-atom exchange (the next largest one is an 
order of magnitude smaller) must be taken into account in any treatment of the 
phase transition of bcc 3 He. 

We gratefully acknowledge the constructive comments and continuous encourage- 
ment of R. C. Richardson. In addition the assistance of W. P. Halperin, D. R. 
Nelson, and M. E. Fisher is much appreciated. 
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Table I. Exchange frequencies and related intermediate numbers for bcc He 
(24.0 cc/mole) . 
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